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(S| ■ Abstract 
(N 

^ \ Projecting a non-Abelian SU{2) vacuum gauge field - a pure gauge con- 

\ structed from the group element U - onto a fixed (electromagnetic) direc- 

' tion in isospace gives rise to a nontrivial magnetic field, with nonvanishing 

magnetic helicity, which coincides with the winding number of U . Although 
0^ , the helicity is not conserved under Maxwell (vacuum) evolution, it retains 

one-half its initial value at infinite time. 



I. INTRODUCTION 



It has been suggested that primordial magnetic fields can develop large correlation 
^ , lengths provided they carry nonvanishing "magnetic helicity" / d r a-b, a quantity known 

^ \ to particle physicists as the Abelian, Euclidean Chern-Simons term. Here a is an Abelian 

gauge potential for the magnetic field b = V x a. If there exists a period of decaying 
turbulence in the early universe, which can occur after a first-order phase transition, a 
magnetic field with nonvanishing helicity could have relaxed to a large-scale configuration, 
which enjoys force-free dynamics (source currents for the magnetic fields proportional to 
the fields themselves) thereby avoiding dissipation 

In this paper we accomplish two things. First we show that configurations of (a, b), 
with quantized helicity, arise from vacuum configurations of a non-Abelian SU (2) vector 
potential. The quantization occurs because (a, b) wind and intertwine in an intricate 
manner. We relate this "winding number" of our Abelian fields to the topological prop- 
erties of the structures in a non-Abelian SU{2) gauge group. Our construction is based 
on lis (5*3), which is relevant to a non-Abelian gauge theory in 3-space. [The construction 
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can also be described in terms of Hopf maps and n3(S'2), which have aheady appeared in 
the hterature |@]; we show how our expressions appear in that formahsm.] Only integer 
winding numbers occur in the mathematical setting of n3(S'3) or n3(S'2); yet we find that 
half-integer windings also lead to interesting Abelian gauge fields. 

Second, we study the time evolution of the magnetic helicity. Since ^ / d^r a • b = 
-2 / dV E ■ b as a consequence of the definition for the electric field, E, time variation 
is determined by the specific physical situation that fixes E ■ b. For plasma physics or 
magnetohydrodynamics the projection of E on b is proportional to the resistivity of the 



medium, and vanishes for infinite conductivity (zero resistivity), see Sec. |^ In that 
approximation, helicity is conserved. However, cosmological electromagnetic fields can 
be expected also to experience evolution in vacuum, where E ■ b 7^ 0. Therefore, it is 
interesting to determine what happens to the magnetic helicity under Maxwell evolution. 
Specifically, we posit that at t = there are magnetic fields with quantized helicity, and 
no electric fields. With this initial data, we solve the Maxwell equations and find that 
at if: = 00 the helicity is precisely half its initial value. The calculation is carried out 
explicitly for two interesting cases (integer and half-integer quanta) and then a general 
argument is given, which requires a regularity hypothesis. 

In Section |T|, we describe how to construct from non- Abelian vacuum fields Abelian 
gauge fields with quantized helicity. In Section we study time evolution. Conclud- 
ing remarks comprise Section where we also speculate on the applicability of these 
mathematical considerations to the SU{2) x f/(l) "standard model", before and after its 
electroweak phase transition. 

II. HELICITY OF GAUGE FIELDS 
A. Chern-Simons structures 

The Chern-Simons number of a non-Abelian gauge potential is given by 

CS{A) = dhe^=\Atd,Al + \r'^A^A)Al) 

= J d^'re'^HiiAd.Ak + lA,[Aj,Ak]) 

= -±^JtT{AdA+lA') (2.1) 

In the second equality, the Ai are elements of the Lie algebra with generators T°- 

A, = A'^T^ , [T^ T"] = UcT'' , tr T"T^ = -5„,/2 

while form notation is used in the third equality: A = A^ da;*. The normalization factor 
is chosen for later convenience, and is also maintained in the Abelian limit, where 



only the bilinear part of (O) survives. Thus our Chern-Simons quantity is times the 



magnetic helicity, which entered physics in the work of Woltier [Q and was subsequently 
further elaborated by many people, including Moffatt, Berger and Field, as well as Arnold 
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and Khesin |^ . Particle physicists have considered both Abehan and non-Abehan Chern- 
Simons terms on 3-dimensional Minkowski space in studies of (2 + l)-dimensional gauge 
theories 0. 

To show how Abehan fields with nonvanishing and quantized helicity arise from non- 
Abelian vacuum configurations, we recall first that under a gauge transformation 

= U-'^AU + U-^dU (2.2) 

where U is an element of the gauge group, the Chern- Simons term transforms as 

C^(A^) = CS{A) + d(tr dU U-'A) + tT{U~'dUf (2.3) 

With sufficiently regular A and U, the second term on the right side integrates to zero 
ip. However, the last term, also a total derivative - although that is not apparent from 
its formula (but see below) - depends only on the group element, and is not damped by 
any fall-off of A. Indeed its value is the winding number of U, which effects a mapping 
from 3-space into the group. 

W{U) = ^ / dr{U-'dUf (2.4) 

W{U) is an integer when U is nonsingular for finite r and tends to ±/ at infinity, for 
then 3-space can be compactified to the 3-sphere, and we assume that the gauge group 
contains SU{2). Thus we are relying on Yi2,{SU{2)) = n3(S'3) = integers. It follows 
therefore that the Chern-Simons number of a non-Abelian vacuum gauge potential, that 
is, one which is a pure gauge 

A = U-^dU (2.5) 

is the winding number of U. 

Let us now specialize to SU{2). Then U involves the Pauli matrices = a'^/2i. 

U = e""^" 

= cos //2 - ia^cD" sin //2 (2.6) 
and the winding number is explicitly seen to involve a total derivative ^ 

WiU) = —^j d(ea6ccD" du' dcD^ (/ - sin /)) (2.7) 
where cD" is a unit vector in isospace: uj"" = cD"/. 



B. Constructing an Abelian gauge field 

Consider next the Abelian gauge potential 1-form a, constructed from a non-Abelian 
vacuum configuration U~^dU by projecting the latter on a fixed direction in isospace, 
specified by a constant unit vector 
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a = iii rfa^U-^AU (2.8a) 

In components, this reads 

ai = fi^A ^ (2.8b) 
where ^ ^ is the vacuum (pure gauge) non-Abehan potential 

U-%U = A^^ (2.9) 

which satisfies 

e'^%At = -le'^'eabcA'jAl (2.10) 
and carries the winding number 

W{U) = / d're^^'ea.cAtA'^Al (2.11) 

Note that a is not an Abehan pure gauge. 

We now show that the magnetic hehcity of a coincides with W{U). The magnetic 
hehcity, with our normahzation, is given by 

= Y^^"^' / dV e'^^'A Id, A I 

l^n^n"! d?rAte''h,,dA''^Ai (2.12a) 



327r 

where ( |2.8| ) and ( |2.1U| ) have been used. But 

e'^^A 1A ]A t = leacde'^'^ea'b'c'A fA J'^ t 

Thus 

H{a) = / d\e'^\atcAtA',Al = W{U) (2.12b) 

and we conclude that the magnetic helicity is quantized by the winding number of the 
non-Abelian vacuum configuration. 

The form of the Abelian potential a may be given explicitly. Parameterizing the unit 
vector tD" as 

tD" = (sin 6 cos $, sin 6 sin $, cos 9) (2.13) 
and taking n"- to point in the third (electromagnetic) direction, we find 

a = cos e d/ - (sin /) sin 6 dG - (1 - cos /) sin^ 6 d-l* (2.14a) 
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An alternate formula presents ( |2.14a| ) in "Clebsch" form, which will be useful later. (The 
Clebsch form for an arbitrary 3- vector V is V = V7 + aV/9 where 7, a, (5 are three scalar 
functions.) 



a = d(-2$) + 2(1 - (sin^ //2) sin^ G) d($ + tan-^[(tan //2) cos0] 
The magnetic field is determined by the 2-form obtained from (|2.14a|) 



(2.14b) 



da = sin 9 



;i - cos /) d/ de - 2(1 - cos /) cos 9 dB d$ + (sin /) sin 9 d$ d/ (2. 15a) 



or from the Clebsch expression (|2.14b|) 

da = -2 d((sin2 //2) sin^ 9) d($ + tan-^[(tan f/2) cos 9]^ 
Finally the magnetic helicity becomes, according to ( |2.14a|) and ( p. 15a] ) 

H{a) = /(I - cos /) sin 9 df d9 d* 

or from ( p.l4b|) and ( p.l5b|) 

H{a) = J d$d((sinV/2)sin2 9) d(tan-i[(tan //2) cos9 



(2.15b) 



(2.16a) 



(2.16b) 



In the Clebsch parameterization, the magnetic helicity is seen to involve integration of a 
total derivative and is therefore given by a surface integral. 

The explicit Clebsch expressions for a and da demonstrate that one may find two 
magnetic surfaces, Sn (n = 1, 2), which satisfy dadS'„ = b • VS'„ = 

5*1 = (sin //2) sin 9 = c 

52 = $ + tan"^ [(tan//2) cos9j = (po 

c, 00 constants (2-17) 

The intersection of these surfaces forms magnetic lines, that is, integral curves of b that 
solve the dynamical system 



dr(T) 

dr 



b r(r) 



(2.18) 



where r is an evolution parameter for the dynamical system. Evidently, for our configu- 
ration this problem in integrable, leading to curves given by 



COS//2 = Vl - c2 cos(<l> - 00 ) 

c 



sin 9 



sin ($ - 0o) + cos2($ - 



(2.19) 
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C. Hopf Mapping 

It is known that the Abehan Chern-Simons term can be related to the degree of the 
Hopf map, which is quantized according to Ha (5*2) 0. The construction of the relevant 
vector potential proceeds in the following manner. Consider a complex spinor u = 
with u*u = + |M2p = 1. Then = u*a°'u is a unit vector and the Hopf curvature 
is defined by 

V = -^e'^''eabcN''d,N''dkN'' (2.20) 
Because U is divergence-free, it can be written as the curl of a potential Oj given by 



—iu*diU. By comparison with (|2.14b|) we find that our potential arises when u is chosen 

as 



_ -(sinV/2)sin2ee2^*'^'^"'[('^"^/2)(cose)] 
(sin//2)sinee-2** 

and the Hopf index coincides with the helicity of the here-constructed (a, b). 

Evidently, the present expressions are awkward when compared to those based on 
n3(5'3), presumably because the latter construction is directly related to gauge fields. 



D. Explicit Expressions 

Henceforth we work with explicit expressions obtained by identifying uo"" in (|2.6| ) with 
the radial unit vector and taking / to depend only on r. This further means that G in 
$ in (|2.13|) are identified with polar and azimuthal angles of spherical coordinates 9 and 
(/>, with ranges < 6' < vr and < < 27r respectively. (A simple generalization would 
be to identify B and $ with integer multiples of 9 and 0.) The magnetic helicity ( p. 161 ) 
then becomes 

H = -^ Hdr^if - sin/) = -^(/ - sin/) (2.22) 
2tt Jo dr 2tt r=°^ 

where we have taken /(O) to vanish. Observe that when sin/(cxD) is nonvanishing H 
is an irrational/transcendental number. When f{oo) is an even integer multiple of vr, 
sin/(oo) vanishes and H is an integer. But note that also an odd integer multiple of tt 
for /(oo) leads to vanishing sin f{oo), and a half-integer value for H. Therefore it appears 
to us that configurations with half- integer magnetic helicity are also singled out. The 
half-integer fields share further good properties with the integer-valued configurations: 
we shall show below that only for these two the winding number does not change under 
the gauge transformation to the Coulomb gauge. 

When reference is made back to the non-Abelian vacuum configuration, which deter- 
mines a, we see from (p.6|) that for the integer windings U > ± J, as expected. For 

' r— »oo 

the half-integer ones, a "hedgehog" asymptote is attained: U > ±i(T ■ r. 
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The vector potential a and the magnetic field b = V x a are neatly described by 
spherical components, which are 0-independent 

a^ = (cos^)/' (2.23a) 

ae = -(sine)- sin f (2.23b) 

r 

= -(sin^)^(l - cos/) (2.23c) 



br = -2(cos^)-^^(l - cos/) (2.24a) 
f 

be = (sine)— sin f (2.24b) 

r 

f 

60 = (sin0)^(l -cos/) (2.24c) 

(The dash denotes differentiation with respect to r.) The Clebsch representation ( p.l4b|) , 
( ^.15b| ) and ( |2.16b| ) gives a clear picture of the helicity. From ( ^.16b| ) we have 



Svr^ J Stt^ J r sine ^ ad) ' 



Since b^ is 0-independent, the integration is trivial, leaving 

(2.25b) 



H = -— / d^ / rdrb, 
A-n Jo Jo 



=2tt 



The integral is over the positive-x {x, z) half-plane and b^ is the toroidal magnetic field, 
perpendicular to that plane. So H measures the fiux of the toroidal magnetic field 
through a half-plane. 

Note that the functions occurring in our Clebsch parameterization are multivalued, 
owing to the presence of the naked azimuthal angle (p and the tan~^. This is as it must be, 
because the helicity is nonvanishing Indeed the multivaluedness of is responsible 
for the nonzero value of the "surface" integral (|2.25| ), which reproduces the helicity. 



When using a in a description of electromagnetic fields, we are effectively in the = 0, 
Weyl gauge. But also the vector potential must be transverse since there are no sources. 
However, the potential in (|2.23|) is not transverse for general /. Transversality may be 



achieved in one of two ways. We may perform a gauge transformation, transforming a 
to a transverse expression a-^. Alternatively, we may choose the function / so that a is 
transverse. We discuss these two possibilities in turn. 

The transversality condition affects only the poloidal components of a (a^, and ae), 
because the toroidal component (a^,) is ^-independent. By using a gauge function pro- 
portional to cos^ times a function of r, one can choose that function so that the resulting, 
gauge equivalent potential is transverse. The components of the transverse potential then 
read 
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= (cos^)F 
aJ = -{sme)^{r^F)' 

a J = -(sin6')-(l - cos/) 



Fir) = 



f^if-sinf)-^J^dr'r'U--sinf) (2.26) 

Of course the magnetic field remains unchanged, but we must still check the helicity 
integral: while the integrand is gauge dependent, the integral changes only be a surface 
term. It may be that the above gauge transformation contributes from the surface. 
Indeed this happens for general /: one finds 

H(a^) = Hia) + — (1 - cos /) sin / (2.27) 

and the helicity of the gauge-equivalent, transverse configuration differs from the original 
expression. However if, and only if, / goes at infinite r to an even or odd multiple of vr, the 
gauge variance vanishes. Thus both integer and half-integer windings are stable against 
this gauge transformation, which renders the potential transverse, but for other, irrational 
windings, the "winding number" loses its meaning because it is gauge dependent |]^. 

The other way to achieve transversality is to impose that condition on the original 
configuration (|2.23| ), thereby determining /. Transversality requires that / satisfy 



r2/" + 2r/'-2sin/ = (2.28) 

Although this equation cannot be solved by elementary functions, it is easily analyzed by 
analogy to a mechanical problem where "time" is Inr/ro- One finds two solutions that 
are regular at the origin, vanishing linearly with r, and tending to ivr in an oscillatory 
manner for large r. Moreover, the scale of r is arbitrary: the solutions are a universal 
function of r/rg and its negative; the positive solution is plotted in Fig. |I|. Evidently this 
transverse potential necessarily corresponds to half-integer winding |T0|. 



III. TIME EVOLUTION 

According to Maxwell's vacuum equations the transverse vector potential satisfies the 
wave equation, which we integrate subject to the initial condition that a is given at time 
i = by ( p.26| ), with a definite choice for /, and that is zero. Thus at initial time, 
there is no electric field and only the magnetic field (|2.24|) is present, maintained for t < 
by a steady current j = V x b. 

jr = 2{cose)^{f - sin/)' 

j, = -(sin^^)-(/-sin/)" 

r 



j^=(sin^)(l[r(l-cos/)]')' (3.1) 
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FIG. 1. Profile of regular solution to ( |2.28| ) [solid line); profile of 2tan which pos- 

sesses the same x = and x = oo asymptotes (dashed line). 

and carrying energy £ = ^ J d^r 

S = sin^ f/2 £ sin 9 [(sin^ f/2 + r VV4) + cos 29 (sin^ f/2 - r^f^A 



= r%{^- cos m- cos f+r^n 

6 Jo 

We determine the time asymptote for the heUcity in two cases 

/ = 2 tan~"^ r/ro 
/ = 4tan~^ r/ro 



(3.2) 

(3.3a) 
(3.3b) 



The former, with / > tt, corresponds to half-integer winding; for the latter / > 271, 

r— +00 r— »o o 

one has integer winding, also the energy density is spherically symmetric; see (p.2|) [|10| . 
The corresponding magnetic fields are 



half-integer winding 6^ = ~4(cosi 



1 



^2 _|_ ^2 



be = 4{sm9)- 



= 4(sin6')- 



r^ + r^)'^ 
rro 



(3.4a) 
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integer winding hr = — 16(cos 
bo = 16 (sin 0) 



(r^ _l_ ^2-j2 
(r^ + Tq)^ 



^^ = 32(sin^)p^ (3.4b) 

The time-evolved fields are obtained by standard Fourier transform techniques, and 
the helicity integral is evaluated as function of time. We find for the two cases 

H = -l + + 3^ + 73^) (3.5b) 



It is seen that for t oo, H attains one-half its value at t = 0. 

While one would like to have an understanding how the localization of the helicity 
changes with time, it does not seem possible to pose such a question in a meaningful way, 
because the helicity density, namely, the integrand that defines H, is gauge dependent, 
and without invariant meaning. Indeed, as we have noted earlier, even the integrated 
helicity can be gauge dependent when the fields and the gauge function survive on the 
surfaces bounding the integration region. (This points to an analogy with the energy 
in general relativity, whose density is diffeomorphism-dependent. Only the integrated 
quantity is invariant and a unique value is determined only after asymptotic conditions 
are prescribed. Furthermore, the energy may be presented as a surface integral and so 
also can the helicity, when the Clebsch parameterization is used for the gauge potential.) 

The result that under Maxwell evolution H decreases to half its value at infinite time 
can also be understood from a general argument. In the Weyl-Coulomb gauge, which we 
are using, the vector potential satisfies the wave equation, which is uniquely solved with 
our posited initial conditions by 

a^(t,r) = / -^^cosHe-'''"a^(k) (3.6) 
J (27r)'^ 



Here a^(k) is Fourier transform of the initial data (|2.26| ) and ^a'^(t,r) = — E(t, r) 
vanishes at t = 0. It follows that 

H = -i j cos' kte'^''k%J{k)dr (k) (3.7) 
= hl Wf^^ ^ cos2kt)e^^'k^~aJ{krar{k) (3.8) 



By appealing to the Riemann-Lebesgue lemma, we can argue that the term involving the 
cosine disappears at large t, owing to destructive interference, leaving half the value at 
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t = 0. However, this step is justified provided the rest of integrand is well behaved at 
/c = 0, which in turn depends on the behavior of a-^(r) at large r. A dimensional estimate 
shows that a large-distance decrease of a-^ faster than 1/r is sufficient, which would mean 
that b should decrease faster than 1/r^. But ( 3.4a| ) exhibits a large r behavior for b of 
order 1/r^, modulated by an angular factor. Since our explicit calculation supports the 
general argument, the angular factor evidently provides sufficient large-r damping. 



IV. CONCLUSION 

Our investigation is based on the connection between the winding number of non- 
Abelian gauge group elements and the Chern-Simons number of Abelian gauge fields that 
are obtained from the former by projection. This mathematical fact suggests a physical 
scenario for the SU{2) x f/(l) "standard model". Before its (first-order) phase transition, 
its vacuum could be populated by pure gauge, vacuum configurations U~^dU of the SU (2) 
gauge group, which carry nonvanishing winding numbers. After the phase transition, 
one direction in isospace is identified with electromagnetism, and the projection of the 
vacuum configuration becomes a magnetic field with nonvanishing helicity. It remains 
to be shown whether the above scenario is energetically stable and can be justified on 
physical grounds. 

Another problem deserving further study concerns nontrivial evolution of an initial 
configuration with magnetic helicity. Rather than using the free Maxwell equations, one 
would rely on the magnetohydrodynamical ones, which make use of Ohm's law to express 
E in terms of b and j. Its nonrelativistic form is 

j=r/(E + vxb) (4.1) 

where r] is conductivity and v is the fiuid velocity, taken to be divergenceless in a fiuid 
of constant density. Inserting this in the Maxwell equation 

dh 

— + V X E = (4.2) 

produces an evolution equation for B 

^-Vx(vxb) = -lvxj. (4.3) 

Further, approximating j by V x b, that is, ignoring dE/dt because it is negligible on 
the relevant time scales, converts the above into 

^ - V X (v X b) = -V'b (4.4) 

which can be analyzed either with v prescribed externally or determined self-consistently 
by its Euler equation, with initial b of the form ( p.24| ). An interesting choice for an 
external v could be a transverse form that carries nonvanishing "kinetic helicity" / d^r v • 
(V X v); for example, what one gets by taking v in the form a of ( |2.23| ) with / solving 
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( p.28| ). For zero resistivity (infinite conductivity) the right side of ( [4 .41 ) is absent - there 
is no dissipation. Then E = — v x b and magnetic hehcity is conserved, since E • b 
vanishes. 



Useful conversations with C. Adam, A. Polychronakos, and D.T. Son are acknowledged. 
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